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Abstract. We present exact formulas for the form factors of local operators in the repulsive 
Lieb-Liniger model at finite size. These are essential ingredients for both numerical and 
analytical calculations. From the theory of Algebraic Bethe Ansatz, it is known that the form 
factors of local operators satisfy a particular type of recursive relations. We show that in 
some cases these relations can be used directly to derive practical expressions in terms of the 
determinant of a matrix whose dimension scales linearly with the system size. Our main results 
are determinant formulas for the form factors of the operators ('k^(0))^T'^(0) and T'^(O), for 
arbitrary integer R, where T*, are the usual held operators. From these expressions, we 
also derive the inhnite size limit of the form factors of these local operators in the attractive 
regime. 


1. Introduction 

Exactly solvable models play a unique role in theoretical physics and in particular for many- 
body systems. Integrability makes it possible to have an analytical control of the underlying 
physics which is in general out of reach. This gives a remarkable opportunity of sharp 
theoretical investigations of important aspects of both classical and quantum mechanics, as 
reported in many excellent textbooks [1, 2, 3, 4, 5]. For a long time exactly solvable models 
have been the ideal theoretical laboratory where physicists could deepen their understanding 
of many body physics, develop new mathematical tools and test the limits of approximate 
methods developed to be applied in the study of more complicated systems. However, 
integrability in the past has been often relegated to the realm of formal mathematical physics 
and thought to have little contact with the real physical world. This prejudice has been 
removed during the last decade mainly thanks to the crucial experimental advances in the 
field of ultra cold atoms, which have made it possible to realize physical systems that can be 
described, with a good approximation, by ideal models (see [6,7] and references therein). One 
of the prototypical exactly solvable models is the Lieb-Liniger gas introduced in [8]. This is a 
one-dimensional model of bosons with point-like interaction, closely related to experimentally 
realized systems of confined bosons in one dimension [9, 10, 11]. Many other integrable 
models have been engineered in cold atomic laboratories in the last decade, such as the Ising 
spin chain [12] and one dimensional Fermi gases [13, 14]. 

Although integrable Hamiltonians can be analytically diagonalized, the computation 
of most of the physical and experimentally measurable observables remains in general an 
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extremely hard task, beeause of the diffieulty in obtaining manageable expressions for the 
form faetors of a loeal operator, whieh are defined as its matrix elements between eigenstates 
of the model. 

In this work we foeus on the exaet computation of form factors of local operators in the 
Lieb-Liniger gas with a finite number of particles in a finite system. Despite of more than 
fifty years of intense investigation following the seminal Bethe ansatz solution of the Lieb- 
Liniger model [8], there are still many physically interesting quantities for which no practical 
analytical expression is available. To the date, only the form factors of the fundamental 
bosonic field [15, 16] and of the density [17] are known in simple enough forms (in a 
sense which will be clarified later) to allow the computation of the equilibrium correlation 
function at least numerically [16, 18, 19, 20]. For the calculation of expectation values of all 
other local operators in equilibrium states (which are particular limits of form factors), other 
complementary approaches are usually exploited [21, 22, 23, 24, 25]. 

Recently, the importance of disposing of simple, practical expressions for the form 
factors of local operators has also emerged in the study of non-equilibrium dynamics of 
isolated quantum systems following a quantum quench. Indeed, it has been shown [26] that 
late times after the quench an integrable systems is locally described by a single representative 
Hamiltonian eigenstate which can be constructed by means of a generalized thermodynamic 
Bethe ansatz. This approach has been successfully applied to the Lieb-Liniger model [27, 28], 
to the XXZ spin-chain [29, 30], to transport problems [31], and to the sine-Gordon field 
theory [32]. For finite times the same approach can be used in conjunction with numerical 
techniques [26, 33, 34] to provide expectation values of local observables as a single sum of 
intermediate states, instead of a double one of the direct approach [35, 36, 37, 38]. However, 
the calculation of measurable observables is limited by the (un-)knowledge of the form factors 
even for the quench dynamics. Furthermore, the knowledge of more general form factors 
greatly helps also in the determination of exact formulas for the overlaps between initial state 
and Bethe states in the quench dynamics which nowadays are known only in very few cases 
[35,27,39, 40,41,42, 43]. 

The organization of this paper is as follows. In section 2 we introduce the Lieb-Liniger 
gas and its solution via the (Algebraic) Bethe Ansatz. For the sake of clarity we summarize 
the main results of our work in section 3, while their derivation is presented in section 4. In 
section 5 we consider the attractive Lieb-Liniger gas and we show how our formulas can be 
simplified in the infinite size limit. Our conclusions are presented in section 6, while technical 
aspects of our work are reported in the appendices. 

2. The Lieb-Liniger model and the Algebraic Bethe Ansatz 

The Lieb-Liniger model describes a system of bosons constrained on a one-dimensional 
system of length L with periodic boundary conditions. The Hamiltonian written in the second 
quantization formalism is [8] 

Hll= f dx , 

Jo 


( 1 ) 
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where c is the coupling constant (which we now take to be positive, c > 0). The field operators 
satisfy canonical commutation relations ['k(x), = 6{x — y). 

The eigenstates of the system have a well defined number of particles, which is conserved 
by the Hamiltonian (1). The model was originally solved using the Coordinate Bethe Ansatz 
[8], but the natural framework for the computation of physical quantities like the scalar 
product of states, form factors and correlation functions is given by the so called Algebraic 
Bethe Ansatz (ABA) and quantum inverse scattering method [5, 44, 45, 46]. 

One of the fundamental objects in the ABA method is the monodromy matrix 


T(A) = 


A(A) B{\) 
CiX) DiX) 


( 2 ) 


where A(A), B{X), C{X), D{X) are operators acting on a reference state that we indicate with 
|0). These operators satisfy a set of non-trivial commutation relations encoded in the Yang- 
Baxter equations. These, in turn, involve another fundamental object of ABA, namely the 
i?-matrix, a 4 x 4 matrix that in our case reads 

f fid^X) \ 

9{d,X) 1 

1 g{p,X) 

V f{d,X) ) 

where empty entries of the matrix are defined to be 0 and where 


i?(A,/i) = 


( 3 ) 


f(\ \ X-p + ic ic 

f{X,E) = —7 -, 9 [X,p) = 


X-p ’ X-p 

The Hilbert space is generated by the Bethe states defined as 


N 




with dual states 


N 


(oinc(x). 

i=i 

while the action on the reference state of the operators A(A), D{X) is given by 
A(A)|0) = a(A)|0) , D{X)\0) = d(A)|0). 

In the Lieb-Liniger model the functions a(A), d{X) are 
a(A) = , d{X) = . 


( 4 ) 

( 5 ) 

( 6 ) 

( 7 ) 

( 8 ) 


We further define for later convenience the function 


r(A) 


a(A) 

^ ■ 


( 9 ) 


The parameters are called rapidities and can take arbitrary values. The Hamiltonian 

(1) can be related to the transfer matrix defined to be the trace of the monodromy matrix 

(2) , r(A) = trT(A) . This makes it possible to define a 1-to-l map between the -particle 
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eigenstates of the Hamiltonian (1) and the states of the form (5) eharaeterized by a set of 
rapidities satisfying the Bethe equations 


N 




Afc — \j + ic 


k = l 


Afc — Xj — ic 


( 10 ) 


If a state of the form (5) is eharaeterized by a set of rapidities that satisfy the Bethe equations 
(10) we eall such a state on-shell, otherwise we call it off-shell. On-shell Bethe states have 
well defined momentum and energy given respectively by 

N N 

p({^j}) = E • r({\j}) = j2 y ■ (11) 

i=i i=i 

Finally, it is useful to define rescaled operators 




( 12 ) 


where d{X) is given in (8). 

The norm of on-shell Bethe states (5) is given by Gaudin’s formula [3, 5, 47], which in 
our notations reads 

N N ^ N2 , 




i=i 


i=i 


j<k 


where 


N 

= 6jk[L + J2 K{Xj, Xi) ] - K{X„ Xk) , 

1=1 

2c 


7F(A,/i) = 


(13) 


(14) 


(15) 


(A - /r)2 -f ' 

It is possible to explicitly define the action of the field operators \I/(0), ^l(O) on the Bethe 
states of the form (5), [5]. In particular, one can use the following commutation relations [5] 

[vl/(0), 5(A)] = -tV~cA{X) , [G(A), 4/1(0)] = *^5(A) , (16) 


to derive 


N 


N 


N 


'i'(0)n 5(Afc)|0) i\fc AfcCi(Afc) 5(Am)|0) , 


k=\ 


N 


k=l 

N N 


m = l 
m^k 


(0| n C{XJAkd{Xk) , 


(17) 


(18) 


k=l 


where 


N 


k = l rn=l 

m^k 


N 




m = l 
m^k 


m=l 

m^k 


( 19 ) 
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Note that the connection between the ABA and the Coordinate Bethe Ansatz solutions of 
the Lieb-Liniger model can be seen explicitly by computing the wave function corresponding 
to a Bethe state of the form (5). Indeed, using the action of the operator ^(a:) on Bethe states, 
and standard techniques in ABA [5], it can be seen that the following is valid 

N 

{0|>l-(a:jv).,,»(xi)nS(A,)|0) = 

i=i 


= i-tyfcfexp 


T ^ 


VgSn 


n(i 

j<k 


ic sgn{xk - Xj) 


k=l / ^ ' 'Pk 

Consider now the following class of form factors of local operators 

N M 


A-Pi. — A 




( 20 ) 


( 21 ) 


i=i 


i=i 


Note that the above expression is equal to zero unless N — h = M — k. The form factor 
in (21) can in principle be computed either by repeated action of the fields (17), (18) on 
Bethe states or through integration of wave-functions in the framework of Coordinate Bethe 
Ansatz. However, using these methods one arrives at formal expressions which in general 
involve sums of ~ N\M\ terms. These expressions are too complicated for both numerical 
and analytical calculations, there existing no available procedure to simplify them. 

In general, the derivation of practical formulas for the form factors of local operators 
like those in (21) for arbitrary N is indeed a very difficult task. In this paper we address 
this problem and our main results are summarized in the next section. The form factors 
of local operators in the Lieb-Liniger model have been also studied in the recent papers 
[25, 48,49, 50, 51]. In Ref. [52] techniques similar to those exploited in the present work were 
used to study form factors of local operators in systems solvable by Nested Bethe Ansatz. 


3. Summary of our results 

In this section we summarize the main results of our work. These are represented by the 
determinant formulas (22), (27), (30), (32) for the operators (4/f (0))^'k^(0) and 4/^(0) (with 
R arbitrary integer). 


3.1. Form/actor o/(4/^(0))^\['^(0) 


The first result is a practical expression for the form factor of the operator (4/f (0))^\['^(0). 
Let be two sets of rapidities satisfying the Bethe equations (10) such that 

Pj 7 ^ Afc, Vj, fc = 1,..., N. Then, our result reads 

N N r ^ 

(oinc(ft)(>i''(o))"'i'"(o)nB(A,)io>=(-1)^^ n 

i=i i=i j,k=i 


N N 


>^nn 


N 


A, 




nw 




i=i 


detAT {6jk -f Ujk) 

- Vf) (K+ - K-) ’ 


( 22 ) 
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where B, C are defined in (12), Xjk = Xj — Xk and 


N 




Ujk — 


pm Xj ih ic 

_^ Xm, Xj ih ZC 

i Ill=liPm-Xj) 


[K{X,,Xk)-K{Xp,Xk)K{Xs,Xj)] , 


J={Px- P;.)" - 4(P;, - P^)(Qa - Q,) + 3(Pa - E^)\ 

N N N 

j=i j=i j=i 


(23) 

(24) 

(25) 

(26) 


and analogously for P^, E^, . In the above equation, K{X,p) is given in (15) and the 

parameters and Xp are two arbitrary complex numbers not necessarily in the set 
Note that the form factor (22) does not depend on Xp and As as it will be proved in the 
following section. 

A formula for this form factor was already given by B. Pozsgay in [25], but there it 
was expressed as the sum of N{N — l)/2 determinants of N x N matrices. In the case of 
equal sets of rapidities some simplifications occur and Pozsgay’s formula remarkably leads to 
an expression for the thermodynamic limit of the expectation value of ('k^(O))^4/^(0) (and 
analogously for ('ki(0))'^'k^(0)). However the formulas presented in [25] are not very 
convenient in the case of different sets of rapidities, and equation (22) is more suitable both 
for numerical and analytical calculations (see section 3.4 for further discussions). 


3.2. Form factor of 


Our second result is the form factor of the operator 4'^(0) for arbitrary integer R. Let 

be two sets of rapidities satisfying the Bethe equations (10) and such that pj Xk, 
\/j = 1,..., N,yk = 1,..., N + R. Then, our result reads 

N N+R N+R 

( 0 inc(ft)»«( 0 ) n e(V)|0) = ^55PT#rTT!(-i)"‘''"'’ 11 + 

j=l j=l ^ j,k=l 


N+R N 


N+R 


^ n n A, - n ('Aj '^ r , 

j=i k=i ^ j=i 


detAT+iz fSjk + ujp 


V+ - Lp 


R,p R,p I \'^R,s R,s 




where again B, C are defined in (12), Aj^ = A^ — A^ and 
Y± _ nm=l Em — Xj E ic 

■TT’A^+i? \ \ _ 1 _ ' ’ 

11™=! =t IC 


U^R) _ 

'^jk - 


]\l=l{.Pm-Xj 


K, - y«.i n"-- A,) 

m+j 


[KiX„Xk)-K{Xp,Xk)K{Xs,Xf] 


(27) 


(28) 

(29) 


In the above expression K{X,p) is given in (15) and Xp, Xg are again two arbitrary complex 
parameters not necessarily in the set As before, (27) does not depend on Xp and A*. 
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Note that the form factor of ^(0) was first computed in [15], and a simplified expression 
for it was given in [16]. Using the techniques discussed in the next section, it is not difficult 
to see that the expression presented in [16] is equivalent to the case i? = 1 of equation (27). 


3.3. Equivalent formulas 

Equation (22) can be cast in different equivalent forms. We give two of them in the following, 
which are derived in Appendix E. 

Eirst, suppose that the two Bethe states have different momentum P\ P^. Then (22) 
can be rewritten as 

N N 


ncte)('f*(o))>=(o)nB(A,)io) 


i=i 


N 


i=i 

N N 


\N+1 


ij 


= (-1) 

6c(Pa - P^?) 


N 


j,k=l j = l k=l ' / j — i 


IKu - 


detw (ijk + ) 


J = i /c = l 'X ■ ' J 

where \p is again an arbitrary complex parameter, Xj^ = Xj — Xk and 


V+ - V- 

P P 


(30) 


Ijii) _ 
^jk - 






T - Xj) 


[KiX„X,)-KiXp,X,)], 


(31) 


and where K{X,iJ.), Vj, J, Px are defined respectively in (15), (23), (25), (26). Note that 
Ujl'^ is the same matrix appearing in the form factor of the density operator ff'f (0)\l/(0) as first 
derived in [17]. 

Consider now the case where the two Bethe states have equal momentum Px = Pp. Then 
(22) can be rewritten as 

N N 


nc(A.j)(>p*(0))=»=(0)ne(V)l0) 


J=1 


N 


i=i 

N N 


Px=Pi. 


= (-l) 


N_^_ii{Ex Efj,)^ 


N 


j,k=l .7=1 k=l ^ J .7=1 






2Nc 

detiv (Sjk + Ufj}^ 


1 = 1 k=l 'X . 

where Xp is an arbitrary complex parameter, Xjk = Xj — Xk and 


V+ - V- 

p p 


(32) 


fj(2) _ 
'^jk - 


- Xj) 


[K{X,,Xk)-K{Xp,Xk)] + 


U.+ - V- 
1 1 


_ Xk) 


(33) 


and where 7f(A, p). Ex are given in (15), (26) respectively. 

We stress again that equations (22), (27), (30), (32) are valid only for sets of rapidities 
{pj}, {Aj} satisfying the Bethe equations (10) and such that pj Xk Vj, k. In particular, our 
formulas cannot be used for the computation of expectation values (corresponding to the case 

Mil = {Ai}f=i). 
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3.4. Numerical checks and discussions 

All the formulas presented in this work have been numerically checked against exact 
computations for a small number of particles. We remind that the form factors (21) can 
be computed, for small values of N, by repeated action of the field (17) on Bethe states. 
We exploited this property and we numerically computed the form factors of \['^(0) for 
R = 1, 2, 3,4,5 and for a number of particles up to iV + i? = 9. This procedure gives the 
same result of our formula (27). Formulas (22), (30), (32) for the form factor of ('k^(0))^4^^(0) 
have been checked numerically for a number of particles up to iV = 9 by comparison with the 
formulas presented in [25] and we found perfect agreement between the two. 

The determinant formulas presented in this work are very convenient for numerical 
calculations. In particular, equation (32) was recently used in Ref. [34] for a study of the 
relaxation dynamics of local observables following a quantum quench in the Lieb-Liniger 
model. In this work, within the quench action method [26], the exact time evolution of 
the normal ordered observable : p^(0) : (where p is the density operator) was computed for 
systems with a number of particles up to = 96. 

A possible (straightforward) application of equation (27) is given by the numerical 
computation of the class of form factors (4/^(0))^'k^(0) through a single resolution of the 
identity, both in the repulsive and the attractive case (see section 5). This can be done very 
efficiently for example using the so called ABACUS algorithm [18, 20, 53, 54] 

The formulas presented in section 3 are also suitable for non-trivial analytical 
calculations as those performed in [50, 51]. In these works the form factors of the operators 
4/(0),\I/f(0)4/(0) were computed in the thermodynamic limit, starting from the corresponding 
finite size formulas derived in [15, 16, 17]. 

Finally, we comment on the fact that equations (22), (27), (30), (32) are valid only for 
sets of rapidities {pj}, {A^} with pj '4j,k. Note that even if two sets {pj}, {Aj} 

correspond to different Bethe states, it might be that pj = Xk for some j, k. However, this will 
not happen in general because of the constrains imposed by the Bethe equations. One case 
where this happens is given by two different parity invariant sets of rapidities with an odd 
number of particles (K being a positive 

integer). In this case, even if {pj} {Aj} the value 0 belongs to both sets of rapidities and 
our formulas don’t apply. This however is a well-known issue true for on-shell form factors 
(see e.g. [16, 17]) and it is not a limitation in practice when computing correlation functions 
since one can restrict to even numbers of particles as done in [34]. 

Form factors for Bethe states corresponding to sets {pj}, {Aj} where pj = Xk for some 
j, k cannot be obtained as the limit Xj —)■ pk directly from formulas (22) (27). In particular, 
our formulas cannot be used to compute the expectation value of the operator (\l'^^(0))^'k^(0). 
This is because all the formulas presented in this work are valid for on-shell rapidities, while 
one should start from an off-shell expression in order to obtain a meaningful result also in the 
limit of coinciding rapidities Xj —)• pk, for some j,k. 
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In this section we derive the formulas presented above. For the sake of clarity, some of the 
proofs will be given in the appendices. We begin with some general considerations. 

From (12), (17) and (18) it is clear that the form factors in (21) depend only on the 
rapidities {Xj}jL^, and on be. the values of r(A) evaluated 

at the rapidities. We can thus define the following function 

N M 

j=i i=i 

In the ABA approach r(A) can be considered as a functional parameter. Indeed one could 
expand the form factor (34) using relations (17) and (18) to obtain a formal expression where 
r(A) is kept as a functional variable. The Lieb-Liniger model results for the form factors are 
recovered by making a specific choice of the function r(A). In particular, we can choose r(A) 
in such a way that 

, r{pj) = , (35) 


consistently with (8). However, one can make also another choice for the function r(A) which 
is as follows. First, define the following function of x, depending on the parameters 


N 




k=l 


Xk — X + ic 
Xi 


X 


ic 


(36) 


Then, one can choose the functional parameter r(A) in such a way that 
r(Aj) = ^a.({A A,) = H • 

M _ , • 

r{pj) = ^M{.{Ek}]Ej) = n ~— - ~- 

Ek - dj - 


(37) 

(38) 


Of course, the two definitions (35) and (37), (38) are equivalent if the sets {Aj}, {pj} satisfy 
the Bethe equations (10), but for arbitrary rapidities they are not. 

Suppose we choose (37) and (38). Then the expression for the form factor becomes a 
rational function of the rapidities only (that is, no dependence on the functional parameter 
r(A) remaining). We can then define 


{Afc}) = GNMiidk}, {Afc}, {r{pj)}, {r(Aj)}) 


_ j->h,k 




(39) 


where is give in (34). We henceforth focus on the function We stress again 

that this function is defined for arbitrary values of the rapidities, even though it is physically 
relevant only for those satisfying the Bethe equations (10). 

We are now ready to present the main ingredient in our derivation of the formulas 
presented in section 3. It is given by the following proposition regarding some fundamental 
properties of the form factors of local operators [5, 25, 48, 49, 55, 52]. 
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Proposition 1 Consider the function defined in (39) (with 0<h<N,0<k<M). 
Then the following properties hold 

(» ; (40) 

(ii) consider prn ^ ^hen the asymptotic behavior of as a function of pm is 

given as follows 


ita {\}) = <1 ,, , ,.,, ff: 

/4m—(-oo I (lyc)hj~ /V—> 0 , 


(41) 


(Hi) consider {Aj}) as a function of pm £ {Hj}^=i- Then it is a rational function 

and its only singularities are first order poles at pm = A^, j = 1,... M; 

(iv) the residues of the form factors are given by the following recursive relations 


X 


> J= 

M 




Q^ybrni A/c) 

4V ivi M N 

n fidT Hrn) n /(Afc) Aj) - Yl /(Aj, Afe) Yl fidm, Pj) 

A _ 1 A _ 1 1 1 


j = l 


i=i 




i=i 

j^m 


7 


•h,k 


' AT 1 ?\/f 


. 1421 


Properties 1 — 4 are well known in the theory of Algebraie Bethe Ansatz. For eompleteness, 
we present their derivation in Appendix A. 

Note that the reeursive relation (42), first derived in [55], was also diseussed in the reeent 
papers [48, 49, 52]. In partieular, in [48, 49] its eonneetion with reeursive relations appearing 
in integrable quantum field theory were investigated. 

Our strategy for deriving the formulas presented in the previous seetion is very simple. 
As a first step we show that properties 1 — 4 of Prop. 1 uniquely determine the on-shell form 
faetors (39). As a seeond step, we prove that the determinant formulas of seetion 3 satisfy 
these properties. Note that a similar approaeh was used in Ref. [52] for the study of loeal 
operators in nested Bethe Ansatz systems. 

The first step of our derivation is given by the following proposition. 


Proposition 2 Let h, k be two fixed integers, with h,k > 0, and consider a family of functions 
{TLn,m = TiN,M{.{Hj}f=i) {Aj}^i)} depending on two integers N, M, with N — h = M — k. 
Suppose that TL n,m satisfies properties 1—4 of Prop. 1 for every N > h, M > k. To be more 
precise, suppose that, for the given value ofh, k the following hold 

(•I -.1. {A,}*.i) = ; (43) 


(ii) lim {A,}) ^ 


0 , h = 0 , 

{is/c)hJ^j,f_(j^^({pj}j:^m-, {Aj}) , h > 0 ; 


(44) 


(in) TLN,Mi{Pj}, {^j}) as a function of pm G {pj}j=i is a rational function and its only 
singularities are first order poles at pm = Xj, j = I, ... M; 


(iv) the residues are given by the following recursive relations 
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9^dmi ^k) 


X 


N M M N 

n Em) n Aj) - Yl /(Aj, Afc) Yl fidm, Ej) 

3 = i J = 1 J = 1 J = 1 

jjtm i^k jjtk iT^rn 


X ELn- 1,M-1 ({Ej}jf^ m: {AjIjVfc) • 


(45) 


Then 


n^MiE^U Mr=i) = J^^^mHeAU {^.}f=i) • 


1 M 


■h,k 


\N 


1 M 


(46) 


We present the proof of this proposition in Appendix B. This proof elosely follows the one 
given by Slavnov in [56] for the scalar product formula between on-shell and off-shell Bethe 
states. 

In order to prove the formulas presented in section 3, all we need to do is then to show 
that they indeed satisfy the properties 1 — 4 of Prop. 1. This is done in the following, where 
we separately consider the cases of the form factors of \h^(0) and of (\I/^(0))^(\I/^(0)). 


4.1. Form factor of'^^{0) 


Define the rational function as 


/. /~\R ^ 


1 ) 


j,k=l 


N+R N 


N+R 


j=l k=l ^ ^ j=l \ r j} ^ 


detjv+R (Sjk + 


R,p ''R,pI \^R,s ^R,s 


(47) 


where are given in (28) and (29). Before showing that (47) gives the correct 

expression for the form factors of 'h^(O), hence proving (27), we briefly discuss some 
properties of the function "Hr. 

The matrix 6jk + UjJP in (47) has a similar structure to those appearing in the formulas 
for the form factors of 4/(0) and 4/^(0)4/(0) as computed in [15, 16, 17]. When studying the 
properties of the determinant of such matrices, a set of identities regarding the sum of rational 
functions are important [17, 57]. For completeness we discuss them in Appendix D. 

Using these identities, it is easy to prove for example that TLr does not depend on the 
parameters Ap and A^. This is done in the following way. Define for convenience 


_ 11^=1 Aj) _ T/+ _ T/“ 

lU=i 


(48) 


For k = 2,..., N + R add to the first column of the matrix 6jk + UjjP column k multiplied 
by /eY'\ From identity (D.l) in Appendix D it follows that the first column becomes 
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proportional to ~ ^Rp- Exploiting the multilinearity of the determinant we get 


detjv+R (^6jk + detN+R j 




1 

1 

1 


where 

= j 

r r^{R) K{Xs,Xi) 

1 + 

if A: = 1 , 

otherwise . 


(49) 


(50) 


Now, for k = 2,..., iV + i?, add to eolumn k of matrix column 1 multiplied by 

iK{Xp, Xk). Exploiting again the multilinearity of the determinant we get 


detAT+R — detN+R j , 

where 


(51) 



^jk + 


^(R) 




AO , 


if /c = 1 , 
otherwise . 


(52) 


In the final expression (52) Xp has disappeared: we conclude that the l.h.s. of (49) and thus 
ELr in (47) are independent of the parameter Xp. 

To show that ELr does not depend on we proceed as follows. Eor j = 2, N + R 
add to the first row of the matrix 6jk + U^k^ row i multiplied by Using again (D.l) 

of Appendix D, we see that the first row becomes proportional to We can then 

follow a similar procedure to the one used before to conclude that ELr does not depend on 
either. In a similar fashion, one can prove that the r.h.s. of (22) is independent of both A^ and 


Xp. 

We stress here that with the above procedure we also immediately see that ELr, as a 
function of the parameter pm, does not have poles corresponding to the zeroes of 
and (Vj^ — Vf^f), since these factors are canceled by the determinant in the numerator in the 
r.h.s. of Eq. (47). 

We shall now show that ELr satisfies properties 1 — 4 of Prop. 1 . We begin with property 
1 namely that ELr{^, {AjjjLi) = {—iy/c)^R\. Using definition (47) our goal is to prove that 

« " nf=,(Uy - vAi) detR(A,„ + ' 


{i^/c) 


c2«-l(/(-l)! ; 


JJ (Xjk + ic)- 


j,k=l 


iVt-^R.r>)(VL-VR.s) 


= {-zv~crR\, 


(53) 


A R,s 


R,s/ 


where and are given in (28), (29). This can be seen by induction. 

The case i? = 1 is trivial by direct computation. Supposing now equation (53) is true for 
7? — 1 > 1 we show that it is also true for R. Using simple manipulations and the identities 
discussed in Appendix D, it is not difficult to see that the l.h.s. of (53) as a function of Xr does 
not have poles, and it is bounded at infinity. Thus, as a consequence of Eiouville theorem in 
complex analysis, it is a constant. In order to determine this constant, we compute the limit 
A/j —>■ oo of the l.h.s. of (53). 









Exact formulas for the form factors of local operators in the Lieb-Liniger model 


13 


From (29) we see that detR{6jk + cc 1/A|. as —)■ cxd. Expand now the 

determinant of the matrix djk + UjjP along the last row using Laplaee expansion: 

R-l 

det^ {djk + + (1 + f/S)detR_i(^H) , (54) 

fc=i 


where Aj is the (iV — 1) x (A^ — 1) matrix obtained by removing the last row and eolumn 
j from the matrix 6jk + Using definition (29) we see that oc 1/A^ as A/j — y oo. 
However, for fc = 1,..., i? — 1 the matrix Ak eontains the eolumn {u [^,..., Sinee 

oc l/(Aii;)^, and = (9(1) for j, k,fR'we see that all the terms in (54) are of order 
0(1/Xr)'^ exeept for the last one. So 

detfi (^6jk + Ujk^^ ~ (1 + ul^l)detR-i{6jk + u\k~^^) , for A/? -)■ cx) . (55) 

From the definition (28), (29) the following expansions are given by straightforward 
calculations 


1 + U 


(R) 


RR 


1 R{R-l) 
^ 2 


+ 0{l/X\,) , 


(56) 


2i 


VIr - yR,R = (- 1 ) ^ VAS) , 

'R 


and also 

R 

(Xjk + i 

j,k=l 


Ufd(vL-v, 


ic) 


R,j 


= i-i)^-Hc 


Kp-Vr,^ (Vr,s-Vr,, 


X 


R-l 

n (Ajfc + ic) 

j,k=l 


11^=1 [yR-i,j yR-i,j 


+ 0{Xl) . 


yR-i,p yR-i,pj ^R-iy 

Using (56), (57), (58) and the inductive hypothesis for i? — 1 we finally have 

( ' /~\R ^ ^ 

(Vc) 


a1™oo c2«-l(i? - 1) 


u R _ 1 p n n yti 


j,k=l 


i=i 


X 


detfi (djk + U^k^ 


Kp-yR.p) [VIs-Vr,s 


= (-^^^)^i?! . 


(57) 


(58) 


(59) 


The first property of Prop. 1 is thus proven. To see that the second is true, it is sufficient to 
observe that 


'^Ri.idj}^ {Aj}) oc — 0 , for /i^ CX) . 

dm 


(60) 


Furthermore, property 3 is easily seen to be satisfied using simple manipulations and the 
identities discussed in Appendix D. 

Finally, we address property 4. It is enough to consider the case Rat —)• X^+r because 
Rr is completely symmetric in the sets {pk} and {Afc} . Define in the following M = N + R. 
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In the limit ^ the last row of the matrix 6jk + beeomes (0,0,..., 0,1) so it is 
straightforward to compute 


lini 'Hniidj}, {Aj}) = g{pN, Am) 


X 


'N-l M-1 M-1 N-1 

fidj^dN) /(Aj,Am) f{dN,dj) 

.j=l i=l J=1 i=i 


(ix/c) 


R 




(R-iy. 


M-1 


.s{N-l){R-l) 


X(-1) 


M-lAf-1 ^ M-1 


X 


n +fc) n n n - r, 

j,k=l j=l k=l ^ j=l 

HR) 


:^detM-i ( + U-f. 


(61) 


Vi - (Vi - vj^,j 

where V^j and are defined in (28), (29) for the sets of rapidities {Aj}^7^- 

We have thus shown that the function T-Lji satisfies all the properties of Prop. 1 and, as a 
consequence of Prop. 2, that equation (27) is correct. 


4.2. Form/actor o/(\l/^(0))^\l/^(0) 

We now show that the determinant formula in (22) satisfies properties 1 — 4 of Prop. 1 . 

The form factor of (\I/^(0))^\[^^(0) corresponds to the case (/i, k) = (2,2) of the general 
formula (39). Looking at property 2 of Prop. 1 wee see that in order to prove the validity of 
(22) we also need an expression for the form factor corresponding to the operator 'pf (O)'k^(O), 
which we give now. Let be two sets of rapidities satisfying the Bethe 

equations (10). The form factor of the operator \l/f (0)\I/^(0) is 

N Af+1 _ Af+1 


nc(ft)i''(o)i'"(o)ne(A,)|o> 

i=i i=i 


. r 




c^/c 


X 


AT+l N 

nn 


A,- 

j=i k=i J 




N^l 

niv;-v,- 


JJ (Xjk + ic) 

j,k=l 

det jv +1 ( 6 jk + U, 


( 1 ) 


jk 


i=i 


Kp-Kp [Ks-Vi: 


(62) 


where as usual B, C are defined in (12) and Xjk = Xj — Xk. In the above 1/j, f/j^^ are given 
respectively in (28) and (29) (for i? = 1) while 


T=-[{P,-P,f-iE,-E, 


(63) 


where Px, Ex are defined in (26). 

Formula (62) can be proven using the same strategy of the previous subsection. In 
particular, one needs to show that (62) satisfies properties 1 — 4 of Prop. 1. It is now 
straightforward to check property 1, while one proceeds in the same way of the previous 
subsection to see that properties 3 and 4 are valid. We address property 2. First, note that the 
following expansions are valid for ^ oo 


nm=l Era Xj ± ic 


11^=1 Am - Xj ± ic ' nm=i Am “ Xj ± ic 


nm=l Fm, Xj ± ic 

EN.^xf+1 X -^+ ^(■‘-) ’ 


rsj 


En oo , 


(64) 
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N 


N-l 


{pm ~ PN {Pm ^j') “ 1 “ ^(1) ) PN ^ OO , 


(65) 


m=l 


m=l 


T ^ p% + 0 {pn) , pn ^ oo . 

Using (64), (65) and (66) it is straightforward to compute 


( 66 ) 



(67) 


where I-L 2 is given in (47) (for i? = 2). As we proved in the previous subseetion 



( 68 ) 


so we see that property 2 of Prop 1 is satisfied by the r.h.s. of (62) whieh is thus eorreet. 

We have now all the ingredients to prove formula (22) for the form faetor of 


(\l/^(0))^\l/^(0): one needs to apply the very same proeedure just used to prove formula (62). 


The steps are the same as before and the ealeulations present no diffieulty so we will not 
present them here. 

5. Form factors in the attractive regime 

In this section we specialise the form factors given in section 3 to the attractive Lieb-Liniger 
model (i.e. c < 0) when one of the two Bethe states is the ground state, since in this case 
major simplifications occur and the determinant can be carried out explicitly along the lines 
of [19]. We start by very briefly introducing the attractive Lieb-Liniger model and we refer to 
[2, 19] for all the necessary technical details. 

We consider the Hamiltonian (1) in the case c < 0 and we define for convenience 
c = — c > 0. With this choice, the Hamiltonian (1) describes a system of bosons with 
attractive interaction. The eigenstates are again in 1-to-l correspondence with the sets of 
rapidities satisfying the Bethe equations which now read 



(69) 


In the attractive regime equations (69) admit string solutions corresponding to bound states 
of bosons [19] . Here we focus on the case where all the N particles of the system form a 
unique bound state, namely a A^-string. The rapidities forming a A^-string solution arrange 
themselves in the following way 


A -f * 2 ^^ + 1 — 2j) -t- i5j , 


(70) 


where A is a real parameter (the center of the string), and where 6j are small deviations 
vanishing in the limit where the length L is infinity and the number of particles N is kept 
fixed. The ground state of the system corresponds to a A^-string centered at A = 0. 
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Following [19] we now wish to obtain a simplified expression for the form faetors 
presented in seetion 3 when one of the two Bethe states eorresponds to the ground state of 
the attraetive model in the limit of vanishing deviations. First, we consider the form factor 
of (^^(0))^\I/^(0). Let be two sets of rapidities both satisfying the Bethe 

equations (69), with {Aj} corresponding to the ground state (instead the rapidities 
do not need to form a iV-string and correspond to an arbitrary excited state of the system). 
In a finite system the deviations 5j will be in general non vanishing and one can plug the 
rapidities (70) directly in the equation (22), since singular terms do not appear. Taking then 
the limit } —)■ 0 the determinant expression gets simplified. We report in Appendix C the 
calculations. The final result for the form factor (with |{Aj}) being the ground state) is 

N N 

{0inc(fiy)('i''(0))">p^(0)ne(^j)|0> = ai) 

i=i i=i 


(- 1 ) 


N 


6 


c2 


P: - AP,Q, + 3{E, + -NiN^ - 1) 


NT^{N) 


nL.(A‘^+f(^-i)- 


where F is the gamma function and P^, P^, are defined in (26). 

Consider now the form factor of \I/^(0). Let be two sets of rapidities 

satisfying (69), with corresponding to the ground state while to an arbitrary 

excited state of the system. The form factor (27) can be simplified in the limit {5^} —0. The 
calculations are again reported in Appendix C and the final result for the form factor is 

N N+R 

(0|nc(ft)>i''>(0) n e(A,)|0> = 

i=i i=i 


N + R _ r^(iV + fl) _ 

nLi +i{N + R- if ' 


Apart from the per se interest, the attractive Bose gas became recently a standard 
tool to solve the one dimensional Kardar Paris! Zhang equation [58] by means of replicas 
[59, 60, 61, 62, 63, 64]. The form factors provided here are the starting point for the 
calculation of observables which nowadays are known only in the limit of infinite time 
[65, 66]. 


6. Conclusions 

In this work we presented exact formulas for the form factors of local operators in the Lieb- 
Liniger model. Our main results are: (i) the determinant formula (27) for the form factor of 
4/^(0) with R an arbitrary integer; (ii) two equivalent formulations for the form factors of 
(\l/^(0))^\I/^(0), i.e. one determinant formula (22) and a set of two determinants (30) and (32) 
which are valid respectively for different and equal momenta. Along the lines of [19], we also 
showed how our formulas can be simplified when one is interested in the form factors between 
the ground state of the attractive model and an arbitrary on-shell Bethe state. Our formulas can 
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be useful for both numerical and analytical calculations, in equilibrium and non-equilibrium 
settings. 

From the theoretical point of view, equation (27) is particularly interesting because it 
shows that an entire family of local form factors can be expressed as the determinant of a 
N X N matrix of simple form, similar to those already known for the form factors of ^(0) 
and \I/^(0)\I/(0) [15, 16, 17]. It would be extremely useful to find a similar representation for 
the form factors of the type ('k^(0))^'k^(0) considered in [25]. 

Note that in principle the method applied in this work could be used to derive manageable 
formulas for the general class of form factors (21). However, as it should be clear from our 
derivation in section 4, our approach requires a preliminary guesswork: indeed as a first step 
one tries to guess, for the desired form factor, an expression that satisfies the properties 1 — 4 
of Prop. 1. The second step, namely the rigorous proof that the properties 1 — 4 of Prop. 
1 are indeed satisfied, is straightforward if the initial guess is correct. Of course, the initial 
guesswork is not blindfold and one is guided by the constrains imposed by Prop. 1. 

The starting point for our initial guesswork leading to formulas (22) and (27) has been 
given by the already known determinant formulas for the form factors of 'F (0) and (0) 'F (0) 
as first derived in [15, 16, 17]. Our strategy has been to look for “minimal” modifications 
of such formulas that could still satisfy properties 1 — 4 of Prop. 1 for more general form 
factors. While this approach has turned out to be convenient for the form factors of 4/^(0) 
and (\I/^(0))^'k^(0), we have not yet been able to follow this strategy for the general case of 
Eq. (21). 

The Lieb-Liniger model can be obtained as a scaling limit of the XXZ model [3, 67, 
68, 25]. It is then natural to wonder whether the methods applied in this work could be used 
also for the XXZ chain. In particular, it would be interesting to have a single determinant 
representation for the form factor of (and so giving a simpler representation than the 

one in [54]). These issues will be addressed in forthcoming publications. 
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Appendix A. Proof of Proposition 1 

We now present a detailed proof of Prop. 1. We begin with property 1. To do this, we first 
prove by induction that 

N N 

>l-"(0)ne(V|0> = (-iv^)"iV!]qr(A,)|0) . 

j=i i=i 


(A.l) 
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The case = 1 is trivial using (16). Suppose (A.l) is true for N > 1. Then using (17) we 
have 


N+l 


>I>''+‘(0) n B(A,)|0) = 


i=i 

Af+l 


N+l 


N+l 


-Wc 


5^Afcr(Afc) J]S(A, 


k=l 


m=l 

m^k 


| 0 ) = 


N+l N+l 


/'N+l 


X 


nXAi,A.,.)|0) = (-!x/S)''+' ( ni-lAi) ) (iV + l)!|0), (A.2) 

i=i 


m=l 


k=l 

m^k 


where A^ is given in (19) and where we used the inductive hypothesis and the identity 

N+l Af+l 

Y,Ilf(^k,Xm) = N + l. (A.3) 


k — l m = l 
m^k 


Equation (A.l) is thus proved. Now, for a set of rapidities {Aj} satisfying the Bethe equations 
(10) we have 


N 

i=i 

so (A.l) becomes 


N 

vl/^(0)n^(A)|0) = (-*v/^)^iV!|0). 

i=i 

Analogously one can prove 


N 

{oino+X't’W)" 


i=i 


(iV^)^Ar!(0| . 


(A.4) 


(A.5) 


(A.6) 


Using (A.5), (A.6) property 1 of Prop. 1 is immediately proved. 

We now address property 2. An explicit proof using the Bethe wave functions was given 
in [52] for the cases {h, k) = (0,1) and {h, k) = (1,1) . For the sake of completeness we give 
here a proof of the general case. This is based on the properties of the scalar product between 
two Bethe vectors [5, 56]. It is known from the theory of ABA that this is a rational function 
of the rapidities and that 

N N 

Mm (A.7) 

f-Lm ^CXD 

i=i i=i 

From (A.7) it is immediately seen by induction and using (17) that 

N N+k 

Mm (01 nc(ft)'l'‘(0) n f5(Ai)|0> = 0 ■ (A.8) 

flm^OO J-J- J-J- 

i=i i=i 

We now prove, once again by induction, that 
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N+h N+k 

Mm {o| n c(f,j)(<rt(o))'‘<i>‘(o) n e(A,,)|o> = 

flm 

i=i i=i 

N+h N+k 

= J] C(/x,)(vI/t(0))"-ivI/^(0) n ^(A,)|0) , h>l. (A.9) 

j=i 

j+m 


The case = 1 is immediately seen using (18) and (A.8). Suppose now (A.9) is true for 
h > 1. Using (18) we have 

Af+fe+l N+k 

(01 n c(ft)('i''(o))'“+‘'i'‘(o) n = 

i=i i=i 

N+h+l N+h+1 N+h+1 N+k 

= E n /(/‘•.ft)( 0 | n C(ft)('i''( 0 ))‘>p‘( 0 ) n e(A.)| 0 )- (a.io) 

t=l 0=1 i = l 7 = 1 

o+t j+t 




We take now the limit pm —^ oo using the inductive hypothesis 

Af+h+l 7V+fc 

lim (01 n C(ft)(>pt( 0 ))''+''i''=( 0 ) n e(V)| 0 ) = 

J=1 i=i 

N+h+1 N+k 

= iv^( 0 | n C(,.^)(»'( 0 ))'‘»'‘( 0 ) n e(A 7 )| 0 ) 

7=1 j=l 

j+m 

Af+/i+l Af+ft+l Af+/i+l N+k 

+(,+ifh n n c(,.,)(®t(o))'-'>i-‘(o)ne(A,)io>.(A,ii) 

i=i 


t=l 0=1 

t + m 

O^t 


J=1 

j^m 


Using now (18) we have 

N+h+1 N+h+1 

{iV~c) Y n n C(ft)('I''(0))'-‘'I'‘(0)ni3(Aj)|0> 


N+h+1 


N+k 


7 = 1 0=1 

t+m o#"* 

O^t 


j = l 

j^m 


i=i 


N+h+1 N+k 

(01 n c(ft)(>pt(0))''>i'‘(0)ne(A,)|0>, 

i=i 


(A.12) 


7 = 1 
j+m 


showing that (A.9) is true also for h + 1. The proof of (A.9) by induction is completed and so 
property 2 of Prop. 1 is true. 

Property 3 is known from the theory of ABA [5, 47] so we finally discuss property 4 of 
Prop. 1. The general case for h,k > 0 can be proven once again by induction. The case 
h = k = 0, corresponding to the scalar product of Bethe states, is well known from the 
general theory of ABA [5, 47, 56]. Using the action (17), (18) of the local fields on Bethe 
states one can proceed by induction as done before for property 2. A detailed derivation is 
also given in Appendix B of Ref. [48]. Since it presents no further difficulty, we will not 
report a derivation here and we refer to [48] for the details. 
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We present here a proof of Prop. 2, following [56]. 

Suppose that the funetion 'Hn,m{{E 3 }^=ii {^j}jLi) satisfies properties 1 — 4 of Prop. 1 
for a speeifie value of {h, k) and for every N > h, M > k. We suppose without loss of 
generality that k > h (the proof is analogous in the ease k < h) and define i = k — h. Note 
that by hypothesis we have M = N + i. Consider then quantity 

= . (B.i) 


We prove by induction on N that 

1, {V}f«) = 0, N>h. 


(B.2) 


This is enough to prove Prop. 2. The base of the induction is N = h and it is true since by 
hypothesis EInm satisfies property 1 of Prop. 1. 

Suppose that (B.2) is true for a given N > h. We show that it is also true for iV + 1. To 
see this, consider Tat+i as a function of pi . Thanks to property 3 it is a rational function with 
at most first order poles at A^, j = 1,..., iV + 1 + f. Now, because both and 

J^N+i,N+i+e satisfy property 4, we can compute the residue at pi = Xk 


rAr+l({/ijr}, {Aj}) '^N,N+£{{Pj}{Xj}j^k) 






(B.3) 


where in the last line we used the inductive hypothesis. Thus Tat+i as a function of does 
not have poles. Furthermore, since 'HN+i,N+i+i and both satisfy property 2 of 

Prop. 1, they have the same limit as pi —)• oo, that is 

lim rAr+i({/ij},{Aj}) = 0 . (B.4) 

Summarizing, TAr+i({/ij}, {Aj}) as a function of pi is a rational function with no poles on 
the complex plane and vanishing at oo. Invoking Liouville theorem in complex analysis we 
thus conclude that it is identically 0 and thus Prop. 2 is proved. 


Appendix C. Form factors in the attractive case 

In this appendix we explicitly derive the formulas for the form factors in the attractive regime 
when one of the Bethe states corresponds to a A^-string, as for the ground state. Our derivation 
follows the one presented in Ref. [19], with only few technical complications arising. We 
discuss in detail the form factor of 4/^(0), the case of (4/f(0))^4/^(0) being completely 
analogous. 

Define M = N + R, and consider the set of rapidities {Xjjji-^ arranged as 
Xj = A + i—(^M + 1 — 2j') i6j . 


(C.l) 
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As a first ingredient we need the following expansion for {5^} —>■ 0 

~ + 0(1) , k = j + 1 , 


'’jj+i 


■, sA7 + 0(1). = 

^(1) > k ±1, 
where we defined = Sj — 5k- We define further 


(C.2) 


^ Y[m=l(drn Aj) 


j — 0,1,M + 1 , 


(C.3) 


n"..(A„-A,) ’ 

where Aq and Am+i are given by equation (C.l). Finally, we need the following expansions, 
which can be easily derived 


V+ Vq , j 1 , 

^ j = 2,...,M, 


= 


J = 1,...,M-1, 


V, 


M+l 


j = M, 


(C.4) 

(C.5) 


where V^j is given in (28). 

We now plug the rapidities (C.l) into (27) and take the limit 0. Consider the 

matrix 

W,k = ~ [K{X„ Xk) - K{Xp, Xk)K{X„ A,)] , (C.6) 

^ ||m-l (Am — Xj) 

so that 




M 


n fh - ''id +df)='"d't mw, 


i=i 


jk 


{C.l) 


The matrix Wjk has singular behavior in the limit {5^ } —)■ 0. Since Ap is an arbitrary parameter 


we can choose Ap = 

/ ^ 

01,2 

V2 

.5i,2 

Vi 

<5i,2 

Vi _ 

.5i,2 ^2,3 

0 

0 

0 

0 

V 0 

0 


Vjk 


0 

0 

0 


0 

0 


— Vm-3 _ _ 

<5m-4,M-3 Sm-2,M-1 

Vm-i 
<5m-2,M-1 

0 


-Wi 

-W2 


Vm-1 


Vm-2 


Wm-2 


<5m-2,m-i 

— Vm-2 Wm-1+Vm 


0 \ 
0 


0 

Vm-1 


<5m-2,M-1 

Vm 


>Vm 

Sm-1,M 


<5m-i,m 
Vm-1 


Sm-i,m / 

where we defined Wj = VjK{Xs, Xj). The determinant of this matrix can be computed 
as follows. First, add the first row to the second, then the second to the third, and iterate 
this procedure until the last row. The last row of the resulting matrix has only one entry 
different from 0. Using Laplace expansion on the last row and the fact that the determinant of 
a triangular matrix is the product of its diagonal elements, we finally have 


detMWjk = i-l)^^-^ 


5m-i,m 


'M-l 

.11 AU, 


'M-l 

11 '". 

. J=2 , 


M 


.J=l 


(C.8) 
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Using (C.7) we can plug (C.8) into (27). One then uses (D.l) of appendix Appendix D and 
the following identity, which is easily derived for rapidities satisfying (C.l) 

M 

J] (A^ - A,) = - 1)!(M - j)! . (C.9) 

m=l 

m^j 

Putting everything together, after simple manipulations we get (72). 


Appendix D. Useful formulas 


In this appendix we discuss identities involving sums of rational functions. These are useful 
when studying the properties of the determinant of matrices like those appearing in (22) and 
(27). The prototypical example is 


N+R pj-AT 




i=i 


- a,) 



(D.l) 


where iT(A, p) and are given in (15), (28) respectively. Identity (D.l) is obtained using 
the residues theorem on the complex function 


9s{z) 


{z-Xs- ic){z -Xs+ ^c) (Ar^^ - z) ■ 


(D.2) 


Indeed the function gg (z) has first order poles for z = Xj, j = 1,..., Xn+r and for = Xs±ic 
while it is easy to see that it has vanishing residue at infinity. Using the fact that the sum of 
the residues has to be zero one immediately arrives at identity (D.l). 


Appendix E. Derivation of equivalent formulas 


In this appendix we show how to derive formulas (30) and (32) from (22). First define the 
matrix as 



Sjk + 


K{\p,\k) 

V+-V- ’ 


J J 


m^j 


K(^Xj , A/^) , 


ifj=£, 
otherwise . 


(E.l) 


Consider the matrix Sjk + Ujk appearing in (22). For j = 1,..., N, j i, add to row £ row j 
multiplied hy Qj/Qi, with Qj defined as 


©J = - Vf > (E.2) 

and where is given in (23). With this procedure one obtains 

detN{6,k + U,u) = {V+ - Vg-)deii,M^ , (E.3) 


where the matrix is defined by 



K{\pM) 
V+-V- ’ 

J 3 

^jk T Ujk t 


ifj=f, 
otherwise . 


(E.4) 
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For j = 1,..., N, j i, add now to row j row i multiplied by iK (A^, Xf'Ej, where 

^ Ul=liEm - Xj) 


n™=l (Xm — Xj) 




Using the above procedure we see that det = det nE^, and thus we obtain 

det^iS.k + U,k) = {V+ - Vf)deii,E^ . 


(E.5) 


(E.6) 


Now consider the matrix 6jk + Ujl^ in equation (30). For j = 1,N, j ^ i add to row 
row j multiplied by Qj/Qi, with Qj is defined in (E.2). Using the the identity 

N y-fN , \ \ ^ 

j=l LLm^jy^rn 'Sj) 

which can be derived with the techniques described in Appendix D, one obtains 

detjv(5,fc + - Px)detNM^ , 

where the matrix is defined as 

EPeEE if j = £ , 


(E.7) 


(E.8) 






(E.9) 


5jk + Ujl^ , otherwise. 

Using again simple manipulations and the multilinearity of the determinant one sees that 

det AT = det AT . (E.IO) 


Putting together (E.6), (E.8) and (E.IO) it is now straightforward to see that (30) is indeed 
equivalent to (22). 

We now derive (32) from (22). Erom (E.6) we see that det at= detArE™ for i m, 
so using (E.IO) we can write 

1 ^ ~ 

+ Uit) = (W - . (E. 11 ) 

^ e=i 

The sum of determinants in the r.h.s of (E.l 1) can be rewritten as follows 

N 

Y, det^(Mi,) = (det^i^jk + - det^(5,fc + uj^)) , (E.12) 

e=i 

where Ujd and are given in (31), (33). The validity of Eq. (E.12) follows from 
U-k = U^k’ + ^EiXp, Afc)/(U+ - V-) and from the fact that Njk = iK{Xp, Afc)/(U+ - Vf) 
is a rank 1 matrix. 

Now, from (E.8) we see that detAr(5jA: + oc {Px — P^), and thus its contribution 
vanishes when we consider Bethe states having the same momentum. Putting together this 
observation with (E.l 1) and (E.12) it is then immediate to obtain (32) from (22). 
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